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Abstract 



Disordered solids are known to exhibit quantitative universalities at low temperatures, the most striking of 
which is the ultrasonic attenuation coefficient Q~ l {bj). The established theory of tunneling two state systems 
(TTLS) in its original form (i.e. without extra fitting functions and parameters), is unable to explain this 
universality. While the TTLS model can be modified, particularly by including long range phonon induced 
interactions to explain the universal value of Q _1 , (a) it is not clear that the essential features of the original 
model that has been successful in explaining the experimental data is preserved, and (b) even if it is, it is 
not clear that the postulates of the original model remain necessary. 

The purpose of this study is to derive the universal acoustic absorption and related quantities observed 
in disordered solids by starting from a many-body quantum theory of unspecified amorphous blocks that 
mutually interact through the strain field. 

Based on very generic assumptions and having no adjustable fitting parameters, the frequency and initial 
state averaged macroscopic attenuation (Q~ l ) of a group of interacting disordered blocks is calculated in 
the low temperature regime (T < w) by a novel "trace method" , which then is iterated up-to experimental 
length scales through a real space renormalization group approach. Then using a heuristic second-order 
perturbation argument, the frequency dependence of Q (u>) is found, and combined with the previous 
result to yield the observed universal values in the MHz-Ghz range of frequencies. 

It is concluded that the TTLS postulates are not necessary in order to explain, at least the thermal 
conductivity, velocity shift and sound attenuation of disordered media in the low temperature regime. 



A famous Zen koan asks, "Two hands clap and there is a sound; what is the sound of one hand?". The 
present work suggests, "to logarithmic accuracy, it does not matter". 
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Chapter 1 

Introduction 



The generic term "disordered solid" (or "amorphous solid") is typically used to describe a large variety of 
condensed matter including metallic and dielectric glasses, disordered crystals, polymers, quasi-crystals and 
proteins. Starting from the pioneering work of Zeller and Pohl [1] it has become clear that the thermal, 
acoustic and electromagnetic response of disordered solids are (1) remarkably different than crystals, and 
(2) are remarkably similar among themselves. To put the problem posed by these two statements in proper 
context, it will be necessary to give a little background on the low temperature properties of both disordered 
and non-disordered condensed matter. 

1.1 A Comparison of Crystalline and Disordered Solids 

According to Debye's theory of crystalline insulators, at sufficiently low temperatures the amplitude r of 
atomic motion remains small enough to consider only the first order term in the inter-atomic force Fi — 
2^ijTj. Thus, the low temperature acoustic and thermal properties of all crystalline insulators are determined 
by a Hamiltonian which can be diagonalized into phonon modes, quantized with uniform increments of E = ui. 
In this framework, the standard text-book procedure [5] to obtain crystalline specific heat is to simply 
integrate over these modes to find the average energy density, 



(E)=J2 wD(w)n(u,0)du, 



(1.1) 



where the summation is over phonon polarizations, and n is the number of phonons at a given temperature 
T = 1//3, 

n = R l - (1.2) 

e?" - 1 



and D is the phonon density of states per volume v, at wave number k 



DM = Kf. 



The upper cut-off ujjj of the integral in cqn( 1.1 1 is the energy of a phonon with wavelength A equal to inter- 
atomic distance aQ. Such high end cutoffs are commonly used in condensed matter theories, including the 
present study (cf. Chapter 2), since they conveniently exclude smaller length scales for which the physical 
assumptions no longer hold. This of course, limits the validity of the theory to energy and temperatures 
below ~ he/ ao. 

In the continuum limit, A 3> ao, the linear dispersion relation 

u>i,t = kcij (1.4) 

holds, where q, t is the transverse t or longitudinal I speed of sound in the limit u>it — > 0. At low temperatures 



T <C Wd, eqn(l.l) can be easily integrated to yield the famous T 3 temperature dependence for the specific 
heat. 

d(E) 12^ 4 /T\ 3 

c - = ir = -ir N UJ • (L5) 

Let us now turn to the thermal conductivity, K. Since phonons are the main heat carriers of heat in 
a dielectric crystal, K(T) can be calculated analogously to the kinetic theory of gases. If I and c is the 
frequency averaged mean free path and velocity of a phonon, 

K » -c v cl (1.6) 

Of course, for a plane wave propagating in a perfect infinite crystal whose atoms interact harmonically, I 
and thus K is infinite. However in practice, the phonons are restricted by the finite geometric scale I = L 
of the sample. Thus, if there are no impurities and the temperature is low enough to neglect anharmonic 



interaction between atoms, from eqn(1.5) we roughly get K ~ T 3 . 

At first sight, the Debye theory seems equally applicable to disordered solids, since these materials are 
known to exhibit well-defined phonon modes, as seen from Brillouin scattering experiments below 33 GHz 
[3] and phonon interference experiments below 500 GHz [4]. Thus, at energy scales for which the thermal 
or acoustic wavelength exceeds ap there is no a priori reason for the thermo-acoustic response of a solid, 



disordered or not, to deviate from that of the elastic continuum described above. 

However, experimentally, disordered solids behave nothing like their crystalline counterparts (for two 
comprehensive reviews, see [5J[3]). The specific heat, measurements revealpQ an super-linear temperature 
dependence of specific heat below To, 

C v oc T 1+a (1.7) 

with a ~ 0.1 — 0.3. Remarkably, at 2bmK this term dominates the Debye phonon contribution by a factor 
of about 1000. Above T c the specific heat raises faster before reaching the Debye limit, which manifests 
as a bump in the C v /T 3 vs T curve (a historical convention used to compare the glassy specific heat to a 
"crystalline reference"). To is a more-or-less material dependent temperature, and typically varies between 
1 - 30K. 

According to specific heat measurements, it is clear that something other than phonons is being thermally 
excited in disordered solids. In the established model of amorphous solids, these are assumed to be "tunneling 
two level systems" . In the present work, we acknowledge the presence of these degrees of freedoms but do 
not make any specific assumptions regarding their nature. A more complete description of both models will 
follow below. 

Turning now to the thermal conductivity, measurements reveal a sub quadratic temperature dependence 

DEI 

KocT 2 ~ p . (1.8) 

with j3 ~ 0.05 — 0.2. Again, around a temperature that empirically coincides with the "bump" in specific 
heat, T), this gradually changes to a temperature independent plateau, it is important to note in the present 
context, that even though C v and K are traditionally fit to power laws, as we will see below, logarithmic 
corrections C v ~ T/ InT and K ~ T 2 InT, lead to indistinguishable functional forms below To) 

The universality of disordered solids go as beyond these "exponents" in temperature dependences: In a 



disordered solid, heat is transported by phonons 8 , allowing the use of eqn(1.6l to obtain the ultrasonic 
attenuation coefficient defined in terms of the phonon mean free path, 

Q- 1 = \/(2ir 2 l). (1.9) 

Below T < Tq, with T) again roughly coinciding with that defined above, the mean free path displays 



an amazing degree of universality [71 151 110): In this regime sound waves travel about 150 times their own 
wavelength, regardless of the chemistry or composition of the amorphous matrix. Above To, X/l rapidly 
increases to a (non-universal) constant of the order one. 

Two other universalities are the ratios [5] of transverse and longitudinal sound velocity and phonon 
coupling, 

-«0.6 (1.10) 

Q 

^«0.4 (1.11) 

Xi 

In the present work we will provide some theoretical justification to these ratios, and use their experi- 
mental values 9J as inputs to obtain the value of Q _1 , and related quantities. 

In addition to these quantitative universalities, disordered solids display a series of very unusual acoustic 
non-linearities that are not present in crystals. For example, above a critical sound intensity I c , the acoustic 
absorption is "saturated" [TT] ; i.e. above I c , the solid becomes transparent to sound: 

1 

A second example is the spectacular echo experiments. These are the acoustic analogs of the magnetic 
response of a spin particle in a nuclear magnetic resonance setup, where the amorphous sample is cooled 
down below T < ui. Then two consequent pulses of equal frequency ui separated by time r s ^> t p are applied 
for a duration of t p and 2t p respectively. The two elastic pulses surprisingly induces a "spontaneous" third 
one, precisely r s after the second pulse. 

It is remarkable that, although less pronounced, the universalities hold for metallic glasses as well. 
However, we limit our considerations in this thesis only to amorphous insulators, since the electron-electron 
and electron-phonon interactions add an extra layer of complexity to the problem. Further, we will focus 
only on the universality of linear acoustic properties of disordered solids in the resonance regime, and the 
thermal conductivity below T c . We will not discuss the nonlinear properties, or the acoustic response in the 
relaxation regime. 

1.2 A Brief Review of the Tunneling Two Level System Model 

The theoretical interpretation of the low temperature data on amorphous materials has for 40 years been 
dominated by the phenomenological "tunneling two state system" (TTLS) model [THUS], which has two 



main assumptions; the first is that, within all disordered solids there exists entities that tunnel between two 
metastable states \L) and \R) (conventionally called "left" and "right", although the tunneling coordinate 
is not specified; cf. below). The tunneling entity is described by a reduced double well Hamiltonian, 



H 




TTLS 



Here A is the double-well asymmetry, A is the tunneling matrix element expressed by the standard WKB 
formula 

A ~ Oe~ c 
£ = dy/2M Q V (1.12) 

M Q is the mass of the tunneling entity, £1 is the ground frequency of each (harmonic) well, d is the separation 
between wells, and V is the double well barier height. The TTLS couples to strain e (and of course, the 
thermal phonons) through the coupling constant 7. 

The second fundamental assumption of the TTLS theory is that the two level parameters of an ensemble 
of TTLS are distributed according to the probability density 

P(A,A , 7) .^p^ 

The 1/Ao factor can be obtained by assuming that the parameter £ in the exponent is uniformly dis- 
tributed. 7 is assumed to be the same for all tunneling entities. 

At the cost of introducing a fairly large number of fitting parameters, the TTLS theory gives an attrac- 
tive explanation to the nonlinear acoustic effects [5] and a specific heat linear with temperature c v ~ T and 
quadratic in thermal conductivity K ~ T 2 [T^J H3] reasonably close to experiment below Tq (see however, 



Fig 5.2 and corresponding discussion in Chapter 5). Moreover it gives the temperature and frequency de- 
pendences of the mean free path and sound velocity[2J [T5] (although the latter requires an additional fit 
function with no physical basis). 

The dominant response of a TTLS to a sound wave in the regime T < u occurs through "resonant 
absorption" . A single phonon interacts with a TTLS with matching energy separation at thermal equilibrium. 



The mean free path is found from Fermi's golden rule for transition rates 1/t. Using I = ct, 

\1- 1 {T,lo) = 2tt 2 Q- 1 =27r 2 ^4'tanhf- N ) (1.13) 

pc 2 \T J 

this equation can be plugged in the Kramers-Kronig principal value integral 



i poo 2/ — 1/ l\ 

c(T) - c(0) = Ac = -p.v. / , { J duj' (1.14) 

7T J U} Z - LU U 



to yield a logarithmic velocity shift for kT ~ Tiuj [T^ 

Ac -f 2 P (uj 



c pc 2 



log (|) (1-15) 



Since T < w is experimentally difficult to attain, this equation is usually tested in the regime T > uj with 
respect to an (arbitrary) reference temperature T , 



c(T) pc 2 b V^o 

Therefore, (with one exception |15j ) most of our direct knowledge regarding Q^ 1 (and thus its universality) 
in the resonance regime comes indirectly, through measuring the velocity shift at high temperatures and 
then obtaining Q _1 from the Kramers-Kronig integral that goes the other way around. 

For higher temperatures T > uj, the dominant mechanism for sound absorption is quite different than that 
described above, through a "relaxation" process. This is when a low frequency sound wave modulates the 
TTLS energy E, causing the population ratio n of excited to unexcited TTLS to differ from its equilibrium 
value 

ncxe"^ E . 

It is not difficult to see how acoustic energy is converted to heat this way: The pulse must do work to widen 
the levels of an excited TTLS, which then is "lost" when the TTLS emits a thermal phonon to decrease the 
overpopulated n. In the limit u>T m i n ^C 1 (which corresponds to higher temperatures), the theory predicts 
half the value of the absorption in the resonance regime, again independent of temperature and frequency 

Q- 1 = ^ (1-17) 

2pc z 

where T min is the minimum relaxation time among the TTLS with splitting E = kT. 



1.3 Problems with the Standard TTLS Model 

The first difficulty of the TTLS model is the uncertainty regarding its microscopic nature. Although the 
configuration coordinate is typically referred to as "left" and "right" , neither the entity undergoing the 
tunneling motion, nor the nature of its motion is known. The TTLS may correspond to a single atom 
translating between two relatively low density regions in the solid, or two atoms may be sliding or rotating 
by each other. A recently discovered isotopic effect [T7] proposes that the TTLS is as large as ten or twenty 
atoms. 

Various mechanisms that could give rise to tunneling two levels were suggested for individual amorphous 
solids, such as the mixed phase metallic NbZr alloy domains|18j or the hydroxyl ion impurities in vitreous 
silica [19] . The tunneling impurities are also known to be present in disordered crystals, such as KBr-KCN 
solutions {20] and are extensively studied. However the following questions remain: First of all, even if two 
level excitations are present in (at least some) disordered media, obviously the are not the only kind of 
excitation that is possible. Why then should the thermo-acoustic response of disordered media be entirely 
governed by two level excitations? And furthermore, even if these excitations are present and dominant 
among all, it is not clear why their distribution of parameters P(A, Ao,7) should be so similar in every 
amorphous material. 

Ideally, a theory for disordered condensed matter should start from assumptions that holds equally true 
for a very wide class of materials. Even if the TTLS hypothesis is one that is effectively correct, it seems 
important to base its assumptions on a firm physical basis. 

Secondly, although the standard TTLS model can explain qualitative features such as phonon echoes and 
saturated absorption as well as approximate temperature dependences of quantities such as the specific heat 
and ultrasonic absorption [3 [5], in its simplest form (i.e. without ad- hoc additional fit parameters and fit 
functions) it fails to explain velocity shift and absorption data at low temperatures [2T1 1221 115] . in addition 
to the the "bump" and "plateau" that appears in the specific heat C V (T) and thermal conductivity K(T) 
data at temperatures above T c [TJ [7J . 

Thirdly, in the temperature regime where the TTLS model is supposed to work (T < IK), it fails to 
explain the quantitative and qualitative universalities [23, 23] the most striking one of which is the ultrasonic 



attenuation Q 1 [TO]; it is not at all obvious [71 HM1 12"3"] that the coefficient of eqn(1.13l 



a=^ (1.18) 

should be material independent. In fact, the quantities in the numerator and denominator are independent 



parameters of the TTLS theory, and if we restrict ourselves to dielectric glasses and polymers, while /!/ = pcf 
varies from material to material by nearly three orders of magnitude, Q~ x remains around (3 ± 2) x 10 -4 . 
It seems extremely unlikely that 7 2 and P happen to correlate with fj, to this extent by pure coincidence. 

Finally, the model, in its original form, neglects the fact that, as a result of interaction with the the 
phonon field, the stress of each TTLS's must be coupled[25]. Although the original paper was formulated 
in the language of the TTLS, its outcome is neither sensitive to the number of levels of the interacting 
excitations, nor to the probability distribution function of their parameters. The only requirement is that 
the stress is coupled to the phonon field linearly (c.f. Appendix- A). 

Recently this interaction was incorporated in the TTLS paradigm [26], and was shown that it leads to 
the experimentally observed small universal value of Q^ 1 . However, it seems legitimate to ask whether the 
original features of TTLS, such as those that give rise to saturation of absorption and echoes will survive 
this significant modification. Furthermore, it may be that the original assumptions of the TTLS model may 
be unnecessary after this modification. It is possible that the presence of elastic interactions gives rise to 
the anomalous glassy behavior, regardless of what is interacting. 

In |24[ 123] it was conjectured that if one starts from a generic model in which at short length scales 
there is a contribution to the stress tensor from some anharmonic degrees of freedom, and take into account 
their phonon mediated mutual interaction, one will obtain the significant features of glasses below IK. The 
goal of this project is to quantitatively justify this conjecture, by calculating the frequency and temperature 
dependence of Q~ l and the quantities related to it in the low temperature "resonance" regime T <C w. 
The relaxation regime, non-linear effects or the intermediate temperature (T > T c ) behavior will not be 
considered in this work. 

The layout of the thesis is as follows: In chapter-2 the precise details of the model is defined, and the 
central object of the study, namely, the dimensionless stress-stress correlation, whose thermally-averaged 
imaginary part is the measured ultrasonic absorption Q _1 is introduced. In section 3 a real-space renormal- 
ization calculation of the average of Q~ 1 (a;) over the frequency to and the starting state m (for details of 
the notation see below) is carried out. Here it is shown that this quantity vanishes logarithmically with the 
volume of the system and for experimentally realistic volumes, and has a small value ~ 0.015. In section 
4, on the basis of a heuristic calculation up to second order in the phonon-induced interaction, it is argued 
that the functional form of <5 _1 (w) at T = should be (lnu;) -1 , and that when we combine this result with 
that of section 3, the numerical value of Q^ 1 for experimentally relevant frequencies should be universal up 
to logarithmic accuracy and numerically close to the observed value 3 x 10~ 4 . In section 5 we attempt to 
assess the significance of our calculations. 



Chapter 2 

The Model 



Imagine the disordered solid as being composed of many (statistically identical) cubes of size L. The precise 
value of L need not be specified, so long as it is much smaller than the experimental sample size, yet much 
larger than the (average) inter-atomic distance a. For this system we can define the strain tensor as usual 

where u(f) denotes the displacement relative to some arbitrary reference frame of the matter at point 
r = (x ll x 2l x 3 ). 

The Hamiltonian can be defined as a Taylor expansion up-to terms first order in strain eij 

H = H Q + Y,e l3 f lJ +0(e 2 ) (2.2) 

ij 

where the stress tensor operator T^ is defined by 

f tj = dH/detj (2.3) 

Note that, in general, in a representation in which H is diagonal, T t j will have both diagonal and off-diagonal 
elements. 

We can define the static elasticity modulus X , a fourth order tensor, as 



( o) _ 1 d{T l3 ) 

XijM ~ L 3 Q e .. 



1 / d 2 H 



L 3 \ deijdeki 

e(T) \ %3 . 



(2.4) 
e(T) 



where the derivative must be taken at the thermal equilibrium configuration. For i > a an amorphous 



solid is rotationally invariant. Due to this symmetry, any component of Xlj-ki can ^ e written in terms of two 



independent constants 27J, for which we pick the transverse Xt and longitudinal xi response; 

Xijlki = (Xi ~ 2 Xt)<Mfc! + Xt($ik5ji + SuSjk) (2.5) 

In the approximation of an elastic continuum, these are related to the velocities c; and c t of the corresponding 
longitudinal and transverse sound waves (of wavelength A such that a -C A -c L) by 

Xl,t = pelt (2.6) 

where p is the mass density of the material. As expected, the continuum approximation leads to Q~ l = 0. 



Clearly, we must take into account a Hamiltonian more general than (2.2) to describe the acoustic properties 
of disordered solids. 

2.1 The Stress-Stress Correlation Function 



To go beyond the continuum approximation, we add an arbitrary "non-phonon" term H'(eij) to eqn(2.2), 

H( eij ) = H el ( eij ) + H' '(ey) (2.7) 

Here, the purely elastic contribution H e i, is given by 

H el ( eij ) = const. + J ^ 3 r^xgL e y(^)e«(r) + ^ p< 2 (r) (2.8) 

ijkl i 

Of course, the second term is meaningful if the velocity field u is slowly varying over distances a. The 
"non-phonon" term H'(eij) is completely general; we neither assume that it is small compared to the purely 
clastic contribution H e i, nor do we identify it with any particular excitation or defect. As above, we define 
the "elastic" contribution to the stress tensor T^- by 

T}f ^J2x%e M , (2.9) 

ijkl 



and as we have done in (2.2) and (2.3) we can define the "non-phonon" contribution to the stress tensor by 



H' = ^ + 5> J T^+0(e 2 ) (2.10) 



'.; 



% = dH'/dea (2.11) 



10 



Note that the strain e%j may be due to thermal phonons, as well as experimental probing strains. Since all 
the contribution to Q~ l comes from the non-phonon contributions, for the sake of simplifying our notation 
we will omit the primes in H' and TL from now on, and used the unprimcd symbols Hq and Ty to denote 
non-phonon contributions. 

We will now define the non-phonon linear response function at scale L. A sinusoidal strain field with 
infinitesimal (real) amplitude eij, 

e i:j (r,t) = e tf [e*^-*'- w * ) + e"*^-"')] (2.12) 

will give rise to a stress response, (Tij) (which in general is complex): 

(T tJ )(r,t) = (ly>e^- wt ) + (3i i )*e- < ^-*'- w *), (2.13) 

The complex linear response function Xij:ki(<l,w) is defined in the standard way 

In practice, we will be interested in the A>a limit. We may therefore work in the linear dispersion limit 

X(w,£) « x{u,u/ci, t ) = x(w) (2-15) 

It is not immediately obvious that the non-phonon response function Xij-.ki(q,u) will have the isotropic 



form analogous to (2.5 1, especially at the latter stages of the renormalization where we will be considering 
amorphous cubes with sizes comparable to the wavelength; however, since it is clear that any complications 
associated with this consideration are sensitive at our arbitrary choice of building-block shape, we will assume 
that a more rigorous (q-space) calculation will get rid of them, and thus assume that Xijki{w) will have the 



same isotropic form as (2.5|, thereby defining "longitudinal" and "transverse" response functions X/,t( w ) f° r 
cubes of size L ^> a. 

We may now calculate the absorption time r _1 = cl , in terms of the imaginary part of the linear 
response function. From Fermi's Golden rule, 

Qa'M = ^2] = ^ylmXaH (2.16) 



11 



where the quantity Imx a (ui) is given explicitly, in the representation in which Hq is diagonal, by the formula 
X$w(«) = J ^<m|T^|n)(n|T w |m)<5(^„ - £ m - «) (2.18) 

n 

where \m) and |n) denote exact many-body eigenstates of Hq, with energies E m , E n , and p m is the probability 
that at thermal equilibrium the system initially occupies state m, 

Pm = ^-e~ pBm (2.19) 

where Z is the partition function. 



Our main objective is to calculate (2.171, which depends on many-body energy levels and stress tensor 



matrix elements. Note that the formula is quite general; substituting TTLS assumptions and parameters in 



it directly gives eqn(1.13l 



As an interesting side note, we point out that the TTLS mean free path is independent of frequency at 
zero temperature, which causes the Kramers-Kronig integral 



to diverge logarithmically at zero frequency. This suggests that the actual zero temperature form of Q a (u) is 
a weakly decreasing function of decreasing u>. We will discuss this matter further in Chapter-4, and propose 
that the actual frequency dependence is ~ log - (U/uS), which is the closest form to a constant that prevents 
( |2.20[ ) from diverging] 



2.2 Virtual Phonon Exchange between Blocks 

Let us define a "block" to be all the non-phonon degrees of freedom in a region enclosed by a cube of 
volume L 3 and consider a large collection of bare uncorrelated blocks as described by the Hamiltonian 



(2.10). Since the strain e^ includes the phonon field, the exchange of phonons must give rise to an effective 



coupling between pairs of block stress tensors. The phonon degrees of freedom are harmonic; therefore the 



Strictly speaking, we require Q 1 (ui) = lin^^g-i- 1/log "*" e U/ui. 
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stress-stress coupling should have the generic form 

H^ = f drf d? J2 A mi(r~ ?)T tJ {r)T kl (?). (2.21) 

M Jv 2 ijkl 

The function Aiju(r — r") is calculated in the paper of Joffrin and Levelut[25]. For "large" r — r' . it has the 
form [6 

pcf 2n\r — r\ 6 

Ajjfei = —(Sji - rijni)5ik + ( 1 1 J [— (Sytfjw - SikSji - kiSjk 

+2>{8ijn k ni + SikUjni + Siin k nj + 5 jk nini + djinrnk + SkirHTij) - 15n,-nfcnfen;] (2.22) 

where ft is the unit vector along r — f. The interaction ceases to have the 1/r 3 form for length scales smaller 
than ro, where 

he 
(V(r )> - — = J7 . (2.23) 

Beneath this length-scale (beyond this energy scale) the interaction becomes oscillatory, and this can be 
taken into account by introducing a cutoff energy level Uq to any integral or sum over energies. In the 
present theory, the cutoff energy associated with the "microscopic" blocks is one of the input parameters, 
and will be discussed further in the next section. 



Due to condition (2.23), the virtual phonon wavelengths considered in any stage of this calculation will 
be larger than the block size. Thus, we will assume that r — r' can be replaced with the distance between 
the center points of two blocks Ri — R%, and that the stress 



Tij w / f l3 {r)df 

J V, 



is uniform throughout a single block. Then, the interacting many-body Hamiltonian for a collection of 
blocks can be written as 

N N 

H N = Y,H ( o S) + £ 5>^(i? 8 - a.)1$T%\ (2.24) 

s=l s,s'=l ijkl 

s<s f 



Eqn(2.24) represents the Hamiltonian Hq of the "super block" (of side ~ N X ^L) composed by the N blocks 
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of side L; in principle, we can redefine the energy levels and stress tensor matrix elements for a super 
block and iterate the procedure untill we reach experimentally realistic length scales. At a given stage of 
renormalization, every pair T^T' 5 ' ' comes together with a A which is proportional to L~ 3 (note the same 
factor in the definition of x). Thus it is not difficult to see that the procedure is scale invariant, and the 
renormalization group equation of Q^ 1 might converge to a fixed point with increasing volume. 

Given infinite computational resources, one could start from arbitrary microscopic forms of Hq and 



To, turn-on the interactions (2.211 between blocks, diagonalize the N-body Hamiltonian iTv, and Tj, 



dHjsi /de\- and iterate the procedure to see whether or not the final forms -ff/v & n d Ty at the macroscopic 
scale are independent of the starting H and T . 

Unfortunately, the number of levels and matrix elements grow exponentially with increasing volume. 
Therefore it is not feasible to solve the problem stated in this form. We will instead narrow down our input 
parameters and focus on few output observables that can be deduced analytically. 

2.3 Input Parameters 

We believe that it is one of the strengths of the present work that our results do not rely on adjustable 
parameters, or the existence of other microscopic (unmeasurable) universal ratios to explain the observable 
one [28, 29, 3Ql[3lJ[32] (though cf.[26]). The only two inputs on which our outcome depends sensitively are 
the ratios ci/ct and Xi/Xt (cf. below for details of the notation) both of which are observed experimentally 
to vary little between different amorphous systems (cf. also Appendix). Our third input ro, which is the 
size of a "microscopic amorphous block" (defined below) only enters into our equations logarithmically. 

2.3.1 Microscopic Building-Block Size 

As discussed in the introduction, the temperature dependence of thermal conductivity and specific heat 
changes at a critical temperature T c ~ IK — 30K. At 10K the dominant phonons have wavelengths of the 
order 50A, and in 23J it is argued that this is just the scale at which we get a crossover from "Ising" to 
"Heisenberg" behavior. Essentially, it is this length scale at which the approximations used in obtaining the 



simple R~ form of A, eqn(2.21 1 breaks down. Therefore, we take the microscopic starting size to be about 
ro ~ 50A, which is still much greater than a. We identify the ultraviolet cutoff Uq with T c , and expect the 
universality to break down as the block size becomes comparable to atomic size. 

Uo is the only quantity in the present work that is not directly measurable. However we should emphasize 
that (a) strictly speaking Uq is not a free parameter, since on a priori grounds we can assign an approximate 
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value to it (b) it can be experimentally obtained indirectly, from thermal conductivity and specific heat 
data through T c — hc/ro- Most importantly, (c) our results, if they depend at all, depend on Uq only 
logarithmically. 

2.3.2 Meissner-Berret Ratios 

The second and third input parameters we will be using are the ratio of longitudinal to transverse speed of 
sound q/q and phonon coupling constants Xt/Xh which are known to be universal among materials up to 
a factor of 1.2 (cf. ref[S], Fig.l and 3). Even though we use as inputs the experimentally obtained values, it 
is not difficult to produce first-principle theoretical justifications for either (cf. below). 



Suppose that the inter-atomic interactions (which reduce to (2.8) for small strains) is due to some 
arbitrary inter-atomic (or inter-block) interaction 4>{r). Namely, 

N 1 

Hei{ea) = const. + ^ </>(r ab ) + - ^ pu?(f) (2.25) 

ab i 

According to the virial theorem, the potential energy expectation value of a harmonic degree of freedom is 
equal to that of the kinetic energy, therefore for the purposes of obtaining the ratio ci/ct without loss of 
generality we will drop the latter, and the const, term. 

(H) = ^2(<t>aj,). (2.26) 

ab 

Further, suppose that the relative displacement u a b of two blocks a and b are proportional to the distance 
r ab between them. By definition of e, it follows that. 

^ab,x — ^xy^ab^y yZ.Zt ) 

u ab,x = &xxTab,x [Z.Zo) 

While this is clearly true for length scales for which r 3> a, in an amorphous structure large deviations from 
this may occur locally for r of the same order as a. However these should even out in the thermodynamic 
limit N — > oo. 



The speed of sound is related to the real part of the zero frequency response function according to (2.6) 

N 



)2 



JY 



XOt,l~ -Q^l^^ab)) 
V a<b 



a<b 



d\r ab \ d(p(r ab ) d 2 <p(r ab ) ( dr ab 






de 2 i;j drab dr ab \ de v 
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(2.29) 

e ij= 



where i ^ j and i = j give Xt and xi respectively. The first term of the right hand side must be zero due 



to stability requirements. Thus all we need to do is to substitute (2.271 into the second term. For a purely 
transverse strain 



r a b = \J{r a b.x + e xy r ab , y ) 2 + r\ by + r 2 z 



(2.30) 



which can be differentiated twice and substituted in (2.29). Letting J2 a b ~~* L 3 J r 2 drdQ 



1 f </>» 2 , / 2 2 , n 

Xat = Y3 I — ~ r dr I rirf.dll 



(2.31) 



where £1 is the solid angle. Similarly, for a purely longitudinal strain 



r a b = y(r a b,x + e xx r ab/x ) 2 + r 2 aby + r 2 z 



(2.32) 



Doing the same as above, we get 



Xoi 



L :! 



f(r)j 



dr 



ridn. 



(2.33) 



Then, the ratio for the speeds of sounds only depend on the angular integrals, 



Ct 



1 

71 



(2.34) 



which is 6% larger than the experimental (average) value. 

The second ratio Im^t/Iinx; is not as trivial to obtain, since these ratios come from the non-phonon 
degrees of freedom. However an argument similar to the above can be made if we consider a second order 
perturbation expansion of the ground state in the interaction V , and differentiate it twice with respect to 
strain e. Then the (short range) interaction <j)(r) must be replaced with the square of the (long range) elastic 



coupling (2.21), and the above argument goes through. 
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2.3.3 Many-Body Density of States 

On general grounds, we can assume that the normalized density of states of the interacting and nonintcracting 
system, can be written as a power series, 

f(uj) = Y, ckUjk ( 2 - 35 ) 

k 

with dominating powers much larger than unity. Note that the density of states F(E) of a composite system 
is given by the convolution of the density of states fi(E) of the constituent objects; 

F(£) = (/l*/2*/3 *...)(£) (2.37) 

where 

/oo 
f(E-oj')f(oj')du,' (2.38) 

-oo 

From which it follows the dominating power of the density of states must be proportional to the number of 



particles if the system is extensive in energy. Thus, (2.35) and (2.361 holds quite generally. 



Finally, while the actual many-body density of states we deduce from the specific heat data (cf . Appendix- 
B), 

f{E) = const.e^/^ 172 (2.39) 

is consistent with the form assumed above, as we will see, none of our results will be sensitive to the precise 
choice of c„ . 

2.3.4 Initial State Dependence 

Finally, we will specify the initial state \m) dependence of the response function, which will be the simplest 
possible choice consistent with our general assumptions, namely the "random form" , 



Xi m) (w) = const.0(£ m +u) = {pcl)Q^ l 6{E m + u) (2.40) 
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Important: This form will not be used until we consider the frequency dependence of Q 1 in Chapter 4 
(cf. eqn. 4.9). Any conclusion we reach in Chapter 3 is not sensitive to this assumption! 



It should be carefully noted that the TTLS form of Xa is not a special case of eqn(|2.40|); this may be 
seen by noting that the form of Q^ 1 given by the latter is approximately, 



Q- 1 ( Wl T) = Qo 1 (l-e-"/ T ) 



(2.41) 



which is different from eqn ( 1.13 1, though not qualitatively so. An important point that should be emphasized 



is, if the non-phononic hamiltonian consisted entirely of harmonic oscillators Q 1 would be independent of 



temperature, which is qualitatively very different than eqn(1.13|. Intuitively, the ansatz (2.41) describes a 



model intermediate between a harmonic-oscillator and the TTLS one, but in some sense close to the latter. 
In principle, with the knowledge of the many body density of states (and thus partition function), a more 
accurate functional form for the initial state dependence can be found. It is our hope that a more realistic 
m-dependence does not qualitatively alter the calculations that will follow in the next sections. 
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Chapter 3 

The Universality of the Average 
Attenuation 



We introduce this chapter by discussing a simpler system (cf. (33] [23]) in which a collection of spin- like 
objects couple to strain through a coupling coefficient 7. As a result of this interaction we will get an 
effective "spin-spin" interaction which is roughly of the form g/r 3 , where 

VI 2 
pc z 

with r\ a dimensionless number of order 1. If the single-spin excitation spectrum P is assumed to be 
independent of energy, it follows on dimensional grounds [53] that 

Poc -. 
9 

In this model, the dimensionless attenuation coefficient is simply 

2pc z 

Furthermore if we include all phonon modes (longitudinal and transverse) this number if reduced by a factor 
of 3. Our purpose in this chapter is to generalize this argument [33] to a more generic model. We will (a) not 
necessarily assume "single particle" excitations, and (b) take into account not only different phonon modes, 
but all tensor components, and show that (a) and (b) alone lead to a surprisingly small value of frequency 
and initial state averaged attenuation, (Q }. 
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3.1 Coupling Two Generic Blocks 

The central quantity we will be interested in this sectionis the frequency and initial state averaged ultrasound 
attenuation coefficient, defined as 

u 



i r 

U J\b *-~ Jo 



where Nb is the number of levels of the block, and Uq is an energy level of the order Uq = hc a /L, but the 
precise value is not essential for our purposes below. 

Our strategy will be to evaluate the quantity M = Tr(V 2 ) in the Ho and Hq + V cigcnbasis. Let us 
start by considering two blocks only, labeled by 1 and 2. We will follow Einstein's summation convention 
for tensor indices. 

V = (Ay W T litf T 2i «J 

= A-ijklA-i'j'k'l'Tl t ijTi ! i>j>T2 J klT2,k'l' 

If Ho has eigenvectors {|«o)} M can be evaluated as 

M = '^Ki ki^i'j'k'i'{rno\T 1 ^T 1 ^ r \jio){no\T2^iT2,k'v\riio) (3.1) 

mn 

Remember that A depends on the relative positions of block 1 and 2. 
Any expression of the form 

I = y^^FjEn - E m )(no\Tij\mo)(n \T k i\mo) (3.2) 

a 

can be written in terms of an integral of xo by inserting unity inside the sum, 

1= / ^F(E n - E m )S(E n - E m - u)(no\Tij\mo)(no\Tki\mo)(hj (3.3) 

J n 

F(u))xom,ijki(v)du> (3.4) 



Doing this twice in eqn(3.1 ) the the trace can be written in terms of xo( 



M = V!V 2 y^ / A ljk iA l , j , k ,i,xo,n 1 ,iji>j>(u')xo,n 2 ,kik'i'(u")duj'duj" (3.5) 
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Here v\ and V2 are the volumes of block 1 and 2 respectively and both the sums and integrals are over the 
whole spectrum. A brief reminder of notation: The first subscript means \ is the response of "noninter- 
acting" blocks, the second ones n\ and n-i denote the eigenstate block 1 and 2 initially occupy. The primed 
and unprimcd i, j, fc, I are tensor components. 

Since an amorphous solid is isotropic all 3 8 components of Xiji'j'Xkik'i' can be expressed in terms of two 
independent constants. We chose the longitudinal xi an d transverse Xt elastic coefficients. 

Xijkl = {Xi - %Xt)<>ij8kl + Xtihkdjl + SuSjk) 



This holds true for any initial state (to avoid clutter we drop the initial state subscript till eqn(3.8l). Let us 

define x, 

LO _Xl 
x + 2 = — 

Xt 
so that 

Xijki = {xSijhi + S ik Sji + 8 a 5jk)xt- 



which simplifies the massive sum of eqn(3.5) into 



MjklM'j'k'l'Xiji'j'Xklk'l' — ^-ijkl^-i' j'k'l' (x5ij5i>j> + Su'Sjj' + 5ij'Sji') (3.6) 

x(x§kiSk'i' +5 k k'5w +5 k i'5ik')x 2 t- (3-7) 

The symmetries of Aijki simplifies matters a bit further, 

h-ijkl = h-kjil = Mlkj = A-klij- 

Substituting the experimental values [9] of Hi/ [H arL d Xi/Xt into A^jki, the entire sum can be evaluated in 
terms of x, 



1 r 
J4 



2^-ijklAi'j'k'l'X0,iji'j'X0,klk'l' ~ 122Q ot = KQ 



a number independent of the relative orientation of the blocks. Thus, eqn(3.5) becomes 



M = K V -^ J2 / Qo^ u ')Qo^'Wdu" (3.8) 

W 12 ni n 2 J J 

where W12 is the volume of the sphere with radius equal to the distance between the two blocks. Notice that 
M is precisely proportional to the square of (Q -1 )- Assuming that the blocks are statistically identical, we 
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can write 



M = K Va ^U 2 N 2 (Q^) 2 (3-9) 

V ab 

Now that we have evaluated M in the eigenbasis of H which allowed us to express it in terms of Qq, Nb 
and Uo, we could also repeat the same steps, this time in the eigenbasis of H = Ho + V, and express the 
trace in terms of the above quantities Q, A/& and U , now modified due to the presence of V. 

M = K^U 2 N 2 {Q~ 1 ) 2 (3.10) 

Note that (Q _1 ) is not necessarily the averaged absorption of the superblock 1+2 (or even related to it by 
a simple numerical factor) , because the definition of the latter involves the squared matrix elements of the 
total stress tensor of the superblock, (T^ ' + T i7 - ) and thus contains terms like (n\T^- \m){m\T^- \n) (where 
\m), \ri) now denote eigenstates of H); while such terms were originally (in the absence of V) uncorrelated, it 
is not obvious that they remain uncorrelated after V is taken into account. We shall, however, argue that on 
average those terms are likely to be small compared to terms of the form |(m|T i7 |n)| 2 , because V involves 
all tensor components of T^ 1 ' and T^ 2 ' while the correlation only involves the same component of T*- 1 ^ and 



T^> . If this argument is accepted, we can identify the (Q 1 ) in (3.10) with the physical inverse absorption 
of the superblock. 



TrV is a scalar, and of course its value is independent of which basis we evaluate it. Equating (3.10) to 



(3.9) 



U 2 {Q^) 2 = U 2 {Q- 1 ) 2 . (3.11) 

To relate {Qq 1 ) to the interacting one (Q -1 ), we must find (U)/(Uq). This will be done using the 



following argument: Let us square our interacting Hamiltonian of eqn(2.24) and consider its trace over the 



same manifold as that defined by the integration limits in (3.1) 



Tr(H 2 ) = Tr(i^) + Tr(H V) + Tr(VH Q ) + Tr(F 2 ) (3.12) 

Since pairs of stress tensors are uncorrelated, we can neglect the second and third terms when evaluating 
the expression in the noninter acting basis. Thus, 

Tr(iJ 2 ) - Tr(H 2 ) = Tr^ 2 (3.13) 
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The Ihs is the change in the variance of energy levels. Since the rhs is positive, we can see that the physical 
effect of the interaction is to spread the energy levels. This equation can simply be written in terms of the 



density of states and eqn(3.10l, 



uj 2 f(uj)duj — 
o Jo 



r/» 



w 2 /oH^ = K^U*N2(Q£)(Q£). 

V ab 



(3.14) 



In general, a macroscopic quantum mechanical system has a density of states roughly given by / (w) oc w n , 
where n is a large exponent proportional to the number of constituent particles. However we can afford to 
be more general and assume 



/o(w) = ^2c 0n uj n 

n 



(3.15) 
(3.16) 



where the dominating terms have power much greater than unity (cf. eqn(3.18)). Note that the experi- 
mentally measured temperature dependence of specific heat C v ~ T is consistent with this form for the 
many-body density of states. In Appcndix-B this consistency is shown, and (although not necessary for the 
present argument) the precise form of c„, as well as an approximate closed form for / is derived. For each 
set of coefficients c„ and cq u the normalization conditions require that the two body system has N 2 levels 
whether they interact or not 



P TTTlo+l 

> can — 

n "0 + 1 

n=0 



P 

£< 

n=0 



u 



n+1 



N* 



(3.17) 



Since both integrals in cqn(3.14| have the form 



1=1 uj 2 f(uj)duj = V" c n 
J o „=o 



jjn+3 

n + 3' 



if the terms for which n + 3 w n + 1 dominate the density of states, we can use eqn(3.17l to write 



77™+ ! 



n=0 



n + 3 



2 /\7-2 



(3.18) 
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Figure 3.1: A superblock of unit length consisting of 8 blocks. 



Thus, eqn(3.14) becomes 



-l\/o-l\ M2rT2r/ VaVb 



N£(U 2 U 2 ) = (Q£)(Q£)NtVtK 



i 1 :. 



(3.19) 



Eqn(3.19) and eqn(3.11) are sufficient to solve for (Q 1 ) = (Q a+b ) 



(Q- 1 ) = 



K 



V a Vb 



(Qoa)(Qo b ) < 



-1/2 



(3.20) 



This equation connects the attenuation coefficient of two non-interacting blocks to that of two interacting 
blocks. 

3.2 Coupling N Generic Blocks: Renormalization Group 



We shall use eqn(3.20) to continue adding blocks till we reach experimental length scales. Let start by 



putting two blocks of side rp next to each other, so that we have a super-block with dimensions 2rp x rp x ro 



(Such as A + H in Fig 3.1 1. For this super-block (sb\). 



VqVb 

V 2 , 
u ab / s6i 



16tt 2 /9 
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and 



(Qs b \) 



9K 



1-1/2 



.(Oo^ 1 ) 2 l ^\ 



Next we combine two sfri's to obtain a larger super-block (sfo) with dimensions 2tq x 2r$ x rp (Such as AH 



BC in Fig 3.1). This time 



V a Vb 



v 



ab / sb2 



16tt 2 /9 



and 



(Q sb2 ) 



9K 9K 

-4- 



1-1/2 



(Q- 1 ) 2 16tt 2 16tt 2 
Finally, we extend one more time to obtain a super-block with the same shape as the original (sb^), but 



with dimensions 2r x 2r x 2r (Such as AHBC + FEGD in Fig 3.1) 



V a Vb 



'(' 



16 



ab / sb 3 



16tt 2 /9 



and denoting Q s ^ = Q 2 ^ and Q — Q L l , we relate the attenuation of a cube with sides 2L to that of a 
cube with sides L 



(Qil 



1 



9K 



(Ql 1 ) 2 

1 
(Ql 1 ) 2 



167T 2 167T 2 

1-1/2 

K 



9K ^ 9K 
4^^ + 16- 



16tt 2 



-1/2 



(3.21) 



where Kq w 150. Eqn(3.21) is our central result. It has the very attractive feature that the value of (Q x ) 
is very weakly dependent on vq = r$, and Qq , since K ^> 1/(Qq 1 ). We now consider the effect of iterating 



the step which led to (3.21 ), by combining eight cubes of size 2L to make one of size 4L; for convenience we 



keep the original definition of the low-energy manifold (E m < Uq), though other choices are also possible. 
Since the only point at which Uq actually enters the result is implicitly in the definition of the "average" in 
Q^ 1 the scale invariant nature of the problem implies that all considerations are exactly the same as at the 



first stage, and we simply recover (3.21) with the replacement of L by 2L and 2L by AL. Continuing the 



iteration up to a spatial scale R, we find 



(Q-'iR)) = [(Qo 1 ) 2 + K log 2 (R/r )] 



-1/2 



(3.22) 



where tq is the linear dimension of the starting block. 



Eqn(3.22) predicts a remarkable counter-intuitive low-temperature effect that would be interesting to 
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test experimentally: As v — > oo, the attenuation vanishes logarithmically (see Fig |3.2| ). While we know of 
no reason why this behavior must be unphysical, in practice we would guess that for any finite ultrasound 
wavelength A, R would be replaced by a quantity of the order A. 



<Q(z)> 
0.10 r 



0.08 



0.06 



0.04 



0.02 



10 15 20 25 30 35 



Figure 3.2: The (numerical) renormalization of (Q). Each increment in z corresponds to an 8-fold increase 
in volume, z — 27 corresponds to r ~ lm 

It is well known that the temperature dependence of the thermal conductivity of amorphous materials 
changes at a temperature of the order of IK] in fact, most such materials show a pronounced "plateau" 
extending very roughly between 1 and 30K. At 10K the dominant phonons have wavelengths of the order 
of 50A, and in [23] it is argued that this is just the scale at which we get a crossover from "Ising" to 
"Heisenberg" behavior (formally, at smaller scales the approximations used in obtaining the simple i?~ 3 



form of A, eqn(2.22) breaks down). Thus, we take the "starting" block size, tq, to be ~ 50A(which is still 



comfortably greater than a). Notice that the result (3.22) depends only logarithmically on 7*0, and thus the 
value of Q^ 1 in the experimentally accessible lengthscales will not be particularly sensitive to this choice. 
Thus, for experimentally realistic values of R we find 



(Q^ 1 )- 0.015. (3.23) 

This value is surprisingly small, and more importantly very weakly dependent on the inputs Qq and ro; 
however (Q) is larger than the experimental values measured in the MHz-GHz range. This is likely due to 
the contribution to the average value, of the rapid increase in Q^ 1 at higher frequencies, as manifest in the 
thermal conductivity data around 10K. To obtain the experimentally observed absorption in the MHz-GHz 
range, we will consider the frequency dependence of Q _1 (a;) in the following chapter. 
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Chapter 4 

The Frequency Dependence of the 
Attenuation Coefficient 



4.1 Level Shifts 

The reader should be warned that the argument presented in this section is rather heuristic and unorthodox. 

Suppose that upon turning on the interactions the change in \ is solely due to the change in the two body 

density of states /o(w). In other words, we will be ignoring the changes in the matrix elements of T. The 

two body density of states f{ui) after turning on the interactions V = AT1T2 can be approximately written 

as 

/M = / fo(x)8(x - uj - A(ui))dx 
Jo 

where 

is the second order correction to the energy level E n . From these two equations, 

Q~» _ /H ^ 1 

Qo\u>) ~ /o(w) ~ |l + dA(«)/dw| 



(4.2) 



We could write A in terms of Q, as described in eqn(3.3) 



AM--*,/ / Q^W-Wq-K^D^ (4 . 3) 



where -E^d + E„ 2 q 2 = u>. The unit step function imposes the ultraviolet cutoff such that \E n — E^\ < U . 
In order to take the integral we must know the initial state |ni) and |n 2 ) dependence of the absorption. For 
this, we use the "random" form discussed in subsection (2.3.4), 

Q- 1 (uj) = 9(oj + E n )Q Q - 1 (4.4) 
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where the attenuation coefficient of an input block Q is a constant. This would happen for example, if 
the level distribution and matrix elements were both uniformly distributed for a starting block. Then the 



integral in eqn.(4.3l can be obtained analytically 

A(u)) = K [-u-U-wlog(w/U )] (4.5) 

Thus, upon turning on the interactions the density of states transforms as 

/M = ^ (4.6) 

This of course, can be generalized to a "single shot" calculation where all blocks within volume R 3 contribute. 
The result is to simply replace the factor K by K log 2 (R/L): 



Q~ H = o — (4-7) 

|X Q - 2 log 2 (i?/L)ln([/oM-l| 
rj const. (m(/7 /u/)) if u < w (4.8) 



4.2 Heuristic Considerations 



While cqn 4.7 seems at least qualitatively consistent with the experimental data (see chapter 5), it is not 



even approximately universal since the constant is inversely proportional to Q Q , and even given a cutoff at 



cj i-v Uq does not satisfy (3.22 1 



There are two obvious reasons form (4.7| cannot be taken to be literally. First of all, the denominator 
of (4.6) is the absolute value of a negative number. Physically this means that the perturbation is so strong 
that the levels are crossing. Of course, we know that this cannot happen due to the "no-level-crossing" 
theorem. We will nevertheless assume that when the perturbation is calculated up to infinite order, the 
overall effect of the higher order terms are small in the level structure, and that the density of states is, at 
least in the limit lo — > qualitatively similar to that given by the first order correction. 

The second difficulty is the high-frequency divergence of x, which indicates that the levels are coming too 
close. While this too can obviously not be literally true (due to level repulsion) there may be some physical 
truth in this divergence too, since we know from the thermal conductivity data, that around T c the phonon 
mean free path rapidly increases by few orders of magnitude, and the divergence we are seeing in the first 
order correction might correspond to this rapid rise. We suppose that higher order terms in the perturbation 
expansion will prevent Q~ x from diverging, and therefore place a cut off to the high end by a large value 
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Q Q (which need not be equal to the microscopic Qo defined above, but presumably of the same order of 
magnitude). Thus, we will consider the following qualitative functional form for the attenuation coefficient 

2~V) = 7i a\ nr, V (°) 

Q + Alog(U/u) 



It is possible to support the ansatz (4.9| further: We know from mean free path measurements, that 
Q~ x is not "noticeably dependent" on to. However if we accept this statement literally, the zero frequency 
velocity shift as obtained by the Kramers-Kronig formula 

Ac_ A rg-vw (410) 



c 



diverges logarithmically. To obtain a finite speed of sound, Q must vanish faster with decreasing to. In 
fact, it can be shown that the closest form to a constant with non-divergent speed of sound is Q^ 1 ~ 
lim epsi j o „_> + l/log 1+e (U/uj). 



We shall therefore accept eqn(4.9) to be the frequency dependence of Q 1 (w), without necessarily as- 
suming "its value" A. 

4.3 The Value of Q- l (lMHz) 



Let us now ask what the value of A must be, in order for eqn (4.9 1 be consistent with result (3.23). The 
frequency average can be evaluated analytically, 

1 f Uo did _ e l/(^«o) 

Here, Ei(x) is a special function, defined as 



Ei(x) = / —dt 



which has the asymptotic form 

e 2Ei(-x) 

lim = = e 7 

where e 27 « 3.172 ... is the Euler-Mascheroni constant. Thus, 



. e l/(,AQ ) / e 2 7 



° ^^A—^^Ql)^ 



2!) 



from which one can obtain the value A ~ 350 that produces (Q) = 0.015. Now that we know A, we may use 



eqn(4.9l to find the value of Q in the universal regime. The experimental probing frequencies u> are typically 



of the order of MHz, therefore we find 



Q(lu = 1MHz) = 2.7 x 10" 



which is precisely the "typical" experimental value. The dimensionless inverse mean free path X/l is then 



X/l « 1/200 



We emphasize that our central result, namely eqn(3.21), concerns the value of Q, which is averaged over u> 



and to symmetrically. Therefore if our assumptions on uj and to were interchanged (or if an assumption pair 
suitably "in between" was used), one would still get similar numbers. 

Let us relax some of our assumptions and see how sensitive the value of Q is. For example, if dominating 
power in the density of states was not a large number, but an arbitrary one, it is not difficult to see that 



this introduces an extra factor between 1/3 and 1 in cqn(3.18), which in turn alters the value of (Q) by a 



factor of about 1.7, very much within experimental variability per material. Note also that eqn(4.11 ) is very 
robust to fluctuations in Q . Namely, 

<1. 



dQ 



Since the equation (4.9) was obtained by heuristically, it is also interesting to see how sensitive the value 



of Q 1 is to the precise details of the functional form. For example, consider more general forms that are 
"near-constant", and repeat the above procedure to find A(s). 



Q-\uj) 



[Qo s +A(s)log(£// W )p 



(4.12) 



The dependence of the value of Q to the exponent s and microscopic value Qq (cf. eqn(4.12 )) is displayed 



in table-1, which suggests that the universality is more general than that required by the precise form (4.9) 







s = 0.5 


s = 0.7 


s = 0.9 


Qo = 


= 0.1 


0.0024 


0.0011 


0.0006 


Qo 


= l 


0.0024 


0.0010 


0.0004 


Qo = 


= 10 


0.0024 


0.0010 


0.0003 



Table 4.1: The dependence of the average value of Q l (\MHz) to the details of its functional form 
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Chapter 5 

Predictions 



In the previous two chapters we deduced the ultrasonic absorption at zero temperature from generic assump- 
tions and arguments. Unfortunately, to the best of our knowledge no experiment has directly measured Q^ 1 
in the regime T <C u. Therefore we will compare our results with indirect measurements from which Q^ 1 is 
deduced. Two such quantities are the thermal conductivity K and the ultrasound velocity shift Ac/c. As 
we have done in chapter 4, we will assume that the renormalizcd ^ m (w) is roughly independent of initial 
state \m) (apart from 9(u> + E m )) and hope that a more realistic m dependence will not qualitatively alter 
the results. 

5.1 Temperature Dependence of Ultrasonic Velocity Shift 

In the TTLS model, the absorption coefficient is calculated by averaging over that of many TTLS. As 
discussed in the introduction, in the "resonance regime" this is given by 

Q- 1 = Q^tanh (^) (5.1) 

where the coefficient Q^f is predicted to be independent of frequency. The resonance behavior of a generic 
amorphous block is rather similar to that of a TTLS; since T <C ui the unexcited distribution of occupied 
levels are well separated from the excited distribution. Let us therefore designate p\ and P2 as the probability 
that the amorphous block occupies a level in the former and latter groups respectively. In thermal equilibrium 
we have, as usual 



(5.2) 
(5.3) 



Pi +P2 


= 1 


Pi 
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The transitions are due to absorption of a phonon and stimulated emission, 



dpi 
dt 



-W12P1 + W21P2 



where the transition rates W12 and W21 are given by Fermi's Golden Rule. From these equations it is trivial 
to calculate the phonon lifetime, 



r-^a-e-"/^^- 1 



,-E m /T 

Z(T) 



Thus, taking into account the assumption regarding initial state inpdependence Q^ can be calculated from 
the phonon mean free path 

Q- 1 (w,T) = Q^ 1 (l- e -"/ T ) 



where now, 






Note that for u ^> T, Q 1 is only logarithmically dependent on frequency, and independent of temperature, 



consistent with experiment [151, and qualitatively similar to the TTLS prediction (5.1 1. We can find the 
velocity shift from the Kramers Kronig relation; 



c pc z 



-A 



Jo w 



(5.4) 



If we insert (5.1 1 into (5.4) we obtain for T 3> w 



£-«?-(£ 



(5.5) 



Note that while this form works quite satisfactorily for low temperatures, it significantly departs from 
experimental data for higher temperatures [34] . The standard way this difficulty is resolved is by introducing 
an additional fit function to the TTLS density of states P(w) = P, so that 

P{u) = P(l + auj 2 ) 



where a is a free parameter. The present model does a better job without such additional fit functions. If 
we treat Q^j as a constant the integral (5.4) can be taken analytically. 
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0.0010 



0.0005 



«T 1 ) 




Figure 5.1: Normalized inverse mean free path (left) and velocity shift (right). The present theory 
(solid) is compared against TTLS without the ui 2 term (dashed) and experiment [IS] [31] (dots). While the 
functional forms predicted by both models are qualitatively similar at low temperatures, the TTLS model 
must use an additional fitting function n(ui) = no(l +au 2 ) for the density of states to resolve the discrepancy 
in fitting Ac/c data (right). 



Ac/c = A(Q_72)[-e" /T Ei(-w/T) - e -" /T Ei(w/T)] 



(5.6) 



which is precisely the same result as ( 5.5 ) in the us <C T limit. In the presence of [ln(£/o/w)] , the integration 



must be done numerically, and is shown in Fig 5.1 and compared with the TTLS result (with unmodified 



density of states) as well as experimental data [T5J [M] . 



K(T)(WcnT 1 K" 1 ) 




1x10 
5x10 



T(K) 



Figure 5.2: Temperature Dependence of Thermal Conductivity. The present theory (solid) K <- 
T 2 \nUo/T is compared against TTLS prediction K ~ T 2 and experiment [35] (dots) below the "plateau". 
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5.2 Temperature Dependence of Thermal Conductivity 

As already mentioned in the introduction, despite the anharmonic degrees of freedom, phonons are the main 



heat carriers in disordered solids. Thus, using the kinetic formula (1.6), 



K(T) » cl h (T)c{T)l ph (5.7) 

= const. T 2 Q(uj,T) (5.8) 

where cf, (T) = const. T 2 is the phononic specific heat and c(T) and l p h are the speed of sound and phonon 
mean free path averaged over frequency. Since at temperature T, the frequency distribution function is 
sharply peaked at uj ~ AK we can evaluate the frequency averages ("dominant phonon approximation"), 

K{T) = const.T 2 ln(£/ /r), 



which is similar to TTLS prediction of K = const. T 2 , but fits the experimental data better (see Fig5.2) 
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Chapter 6 

Conclusions 



Our main goal in this work has been motivated by providing an explanation to the universal acoustic 
response of disordered solids at low temperatures, which for us entails a description that does not use ad-hoc 
fit functions and parameters or postulate other (unobservable) universal quantities or entities. 

We believe that we have been successful in doing so, in that we have shown that (a) simply starting from 
arbitrary uncorrelated blocks and coupling them elastically yields a (Q^ 1 ) factor that is universal, in the 
sense that its precise value depends sensitively only to measurable quantities Ct/ci and Xt/xi, both of which 
fluctuate only by a factor of 1.2 across different materials, (b) That given the ansatz derived in chapter 4 for 
the frequency dependence of Q(u>), the experimentally probed Q(lMHz — 1GHz) is close to the measured 
value, (c) Other quantities related to Q _1 , namely the temperature dependences of the mean free path, 
thermal conductivity and velocity shift, are consistent with our generic model. Finally, we note that the 
<Q- 1 > 

0.015 - 



0.010 



0.005 



<Qo" 1 > 



10" 4 0.01 1 100 

Figure 6.1: (Q^ 1 ) vs (Qq 1 ) ■ Notice that the former can never be larger than O(10~ 2 ) 

essential ingredient of the present model is the existence of non-phonon degrees of freedom that couples the 
the strain field linearly. Thus, our arguments should apply equally well to disordered crystals (cf. |32j). with 
a slightly different angular dependence of the coupling coefficient h-ijki- That being said, our conclusion 
does not imply that all disordered solids must have the same Q" 1 : If the attenuation coefficient of the 
"microscopic blocks" are significantly smaller than the canonical value, then so will be the macroscopic one. 



However, a Q l appreciably larger than the canonical value would tell against our hypothesis (see Fig 6.1 ) 
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Appendix A 

Derivation of the Phonon-Induced 
Stress-Stress Coupling 



For the sake of completeness we give the derivation of the phonon induced stress-stress coupling defined in 
section 2 (also, cf. [6] and [25]). Consider two generic "defect modes", with (non-phonon) stress tensors T\ 
and Ti residing at points T\ and r*2, such that they are coupled to the phonon field at eij(fi) and ey(^j) 
linearly. The full Hamiltonian, includes three terms: The defects Hd, the phonons H e i, and the coupling H c 



H — Hd + H e i + H c 



\Pm\ 2 , mw w 



II n 



+ y^ e a p(?i)T\ a p + e a p{?2)t2, a p (A.l) 

afi 

Where Hd corresponds to the internal energy of the defects, Ui a (f) is the a vector- component of the dis- 
placement of the atom at position r, of defect i, the subscript fi indicates the phonon branch and u qfJ _ is 
the Fourier component of the phonon mode with momentum q. Namely, if e is the unit vector along the 
direction of q, 



a /~\f i-^i u qv c <wp c 



In the formal treatment of phonons in a crystal, N stands for the number of atoms in a unit cell. While 
defining such a cell (or for that matter defining phonons) is rather questionable for a glass, we may think 
of N as an (arbitrary) number such that the above equation is meaningful for phonons with wavelength 
satisfying Na 3> A. Let us write H c explicitly from the definition of strain; 



H r 



jL^r^q^^e^) (A.2) 



'//' 



■M> 



The effect of H c is to shift the equilibrium positions of atoms such that phonon modes can be defined with 
respect to this equilibrium. We shall define a new phonon coordinate for each q and fi, 

^qfl ^q^ ^OqfJ, 

in terms of this new equilibrium coordinate uo qfl , which can be found from setting 

/ dH 



\du, 



q)i 



= 0. 



^2 mug^uoqfi = -r=^2iqpe qii , a (T ha pe l f- Sl + T 2 . a[3 e^- S2 ) 



That is 



Thus, 

Then the full hamiltonian can be written in terms of this coordinate, 

*-~v I V 2 ,,,, mwl „ rnui^,, „\ 

H = H * + E £ + -2^i«^i - -^kj 2 

Note that the first two terms in the paranthesis is formally identical to the phonon Hamiltonian, whereas 
the third, V c explicitly depends on relative defect positions. That is, 

Vc = - E m ^ q,i f^ n (Ti^e*** 1 + T^ aSi e^){T Xa ^ + T 2 ^^) 

ZUILU 1\ 

Let us substitute u> 9M = gc m , keeping in mind not to sum momentums with q > 2it/a. Separating the 
position dependent and constant terms V c — const. + V(R) with R = f\ — r*2, we find 



V c = const. - J2 ^ ^l^ A^f^s cos(q.R) 

q/i M 



Of course, constant part is rather uninteresting (Hd could as well be redefined to include these constants, 
so that the defects are "dressed" with phonons). It is the position dependent term that mixes the energy 
levels of the blocks, and correlates their stress matrices. 
For a longitudinal mode /j, = I, 

_ Qa 
&ql,a I I 

kl 
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whereas for transverse modes \x — t\ % t<i, 

&qti,aQa &qt2,aQa U 



and 



/*=fit a |y| 



These can be substituted in V(R) when summing over fi 



V(R) = E q - m ^f 1 . a[j f^ sC os(q.R){4 + \) (A.3) 

E ^^T^f^s cos(q.R) (A.4) 



Letting 

we may evaluate V(i?). Ommiting the angular dependences 

q a q0q-r<ls/q A 

and 

q a qs/q 2 

it can be seen immediately that 

[>2tT / Ft. rfi rji 

t> / 2rf, rf, -ll,apJ-2,-yS 

V = / g Ti ta pT 2n 8 oc — . 



The angular dependence can be taken into account to give Eqn. (2.22) 
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Appendix B 

Many-Body Density of States of 
Disordered Solids 



This appendix connects assumption (2.35 1 to experiment, particularly the temperature dependence of specific 
heat C V (T). We also find the density of states the many-body density state of amorphous solids in closed 
form. 

The amorphous specific heat, as found from experiment is, 

a(/?) = |. 

The average energy is related to the partition function as 

U = -d logZ (B.l) 

Since the specific heat of glasses are linear in temperature, the U must have the form, 

U = aT 2 + b = J, + b 
At zero temperature the block is in its ground state, which we (arbitrarily) chose to be zero. Thus 

6 = 0. 



Then from eqn(B.l ), 



Z = Ae a/I3 (B.2) 

Here the constant A is the degeneracy of the ground state. This is because 

lim Z = lim V e~ fiEn 



(3— foo j3— foo 



:W 






is equal to number of states for which E n is exactly equal to zero. We can now find the density of states 
f(E) from 

/>oo 

Zp = / f{E)e-? E dE 
Jo 



by substituting eqn(B.2) on the lhs and expand both the partition function and the density of states, 



Ae a " 3 = A(l+ J + ...)= / [/ (1) (£) + f (2) {E) + . . ]e-? E dE. 
P Jo 



The first term of the left hand side can be obtained if 



whereas the second can be obtained if 



fW{E)=A5(E), 



f ( - 2 \E)=aA. 



This way we can get f(E) up-to all orders. The result is, 



f(E) = A[8{E) + S(E)} 



where 



s(E) = y: 



E n - l a % 



^ n\(n- 1)!' 



(B.3) 



One may ask why there is a delta function in the density of states. This is because, for e close the level 
spacing, the integral 

^ f(E)dE 



must give A, the number of levels at E = 0, consistent with our choice of Eq and eqn(B.2 1. Let's look at for 



which value of n is the most prominent in cqn(B.3|. We will assume that k 3> 1, since these are the typical 
energies accesible to experiment. 



Sk = log 



c k E k-\ 
k\{k-l)\ 



k\ogc+ (k - l)log£-logfc! -log(fc- 1)! 
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Using Stirling's approximation, \ogkl = k\ogk — k + 1 



log S k = fc(log c + log E) - log E - fc log fc + 2fc - 3 - k log(fc - 1) + log(fc - 1) 



differentiating with respect to k and setting it to zero, we find 



k{k- 1) =cE. 



Thus, at energy E, the term k « ycE 3> 1 dominates the density of states. 



S^ = 2VcE-2 



One can find a closed form for f(E) by using an asyptotic form of the modified bessel function of the first 
kind, I v (z). For z ^> f , we can approximate 



h(z) 



<2lTZ 



since the density of states can be written in terms of this special function 



1(E) = A 



^/i(2v / cS) 



E 



5(E) 



(B.4) 



we can simplify 



„2Vc.E 



■^)^ 1/4 ro 



74^1/2 



(B.5) 



This asymptotic form of the modified bessel function of the first kind is well known. It is derived from a 
saddle point approximation, by writing 



J^ Sk ~ / S k dk = I i 

k 

and expanding lnS^ at its sharp peak up to second order. 



e lns «dk 
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